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Flat Central Density Profile and Constant DM Surface 
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ABSTRACT 

The scalar field dark matter (SFDM) model proposes that galaxies form by condensa- 
tion of a scalar field (SF) very early in the universe forming Bose- Einstein Condensates 
(BEC) drops, i.e., in this model haloes of galaxies are gigantic drops of SF. Here big 
structures form like in the LCDM model, by hierarchy, thus all the predictions of the 
LCDM model at big scales are reproduced by SFDM. This model predicts that all 
galaxies must be very similar and exist for bigger redshifts than in the LCDM model. 
In this work we show that BEC dark matter haloes fit high-resolution rotation curves 
of a sample of thirteen low surface brightness galaxies. We compare our fits to those 
obtained using a Navarro-Frenk- White and Pseudo- Isothermal (PI) profiles and found 
a better agreement with the SFDM and PI profiles. The mean value of the logarithmic 
inner density slopes is a = —0.27 ± 0.18. As a second result we find a natural way 
to define the core radius with the advantage of being model-independent. Using this 
new definition in the BEC density profile we find that the recent observation of the 
constant dark matter central surface density can be reproduced. We conclude that in 
light of the difficulties that the standard model is currently facing the SFDM model 
can be a worthy alternative to keep exploring further. 

Key words: cosmology: dark matter - cosmology: observations-galaxies: fundamen- 
tal parameters 



1 INTRODUCTION 



In the A Cold Dark Matter (ACDM or LCDM) model, 
known as the standard model of cosmology, the formation 
of structures in the universe is through a hierarchical pro- 
cess of growth of structures, meaning that small structures, 
like small haloes of galaxies merge to form bigger ones, like 
galaxy clusters and superclusters haloes. In this picture, the 
universe contains around 96 per cent of an unknown form 
of energy that is usually called dark matter (DM) and dark 
energy. The ACDM model can successfully describe cosmo- 
logical observations such as the large scale distribution of 
galaxies, the temperature variations in the cosmic microwave 
backg round radiation and the recent acceleration of the uni- 
verse (IColesll2005l : lLahav fc Liddle1l2010l : iPeebles et ai1l2009l; 
iGuo et al.ll20llh . 

However, recent observations in far and nearby galax- 
ies have shown that the model faces serious conflicts when 
trying to explain the galaxy format i on at small scale s 
llsee Robles fc Matos 2011 for a revie 3; iFriedmannl |20H| ). 
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For instance, in the ACDM simulations the halos present 
rising densities towards the central region beh aving as 
p ~ r _1 well within 1 kpc (Na varro et al.|[2010l ). On the 
other hand various observations suggest that the rotation 
curv es are more consistent with a constant central den- 
sity l|Kuzio de Narav et al.ll2008i ; Ide" Blokll2010l') . this is mos t 
commonly known as the cusp/core problem (de Blok 2010). 

Studying a wide range of galaxies of different morpholo- 
gies and with magn itudes in the interval —22 ^ Mb —8 
iDonato et al.l |_2009) fit the ir rotation cur ves (RC) using a 
Burkert profile for the DM l |Burkerdll995h and found that 



log(/i /M @ pc 



2.15 ±0.2 



remains approximately constant, where 



Mo = poro 



(1) 



(2) 



with po the central DM densi ty and rp the core radi u s. Sim - 
ilar results where found in iKormendv fc Freeman! l|2004l ): 
ISpano et al.1 (|2008|). Exploring furthe r the constant value of 
Gentile et al.l (|2009l ) found that within ro 



/to 



for the DM, 



the DM central surface density in terms of the mass inside 
it (M <ro ) is < £ > ,dm= M <r . /7rro 
gravitational acceleration due to DM felt by a test particle 



72±^M @ pc~'- i ,the 
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at the radius ro was found to be 

g D M(r ) = Gvr < E > ,dm= 3.2±iJ x lO^cms" 2 , (3) 
additionally they reported the acceleration due to the lumi- 
nous matter at ro to be gbar( r o) = 5.7^2 t x W~ 10 cms~ 2 . 

In the ACDM model the galaxies have evolved through 
numerous mergers and grown in different environments, the 
star formation and basic properties of the galaxies are not 
expected to be a common factor among them. Therefore 
explaining both the constancy of no and the core in the 
central regions of galaxies seems very unlikely in this model. 

These problems can be used to test alternative DM 
models. There are other models that do not include DM 
but instead modify the Newtonian force law F= ma, one 
of them is call e d MOND and was proposed by Milgrom 
(|Milgroml |20ld : ISandersI 2009), in this modification New- 
ton's law reads F = ma^i(a/ao), where the fixed acceleration 
scale ao divides the Newtonian and MONDian regimes, for 
x « 1 ft(x) = x we have the MONDian regime and for 
x » 1 we recover the usual Newtonian acceleration, the 
value of a ~ 1.2 x 10~ 10 m/s 2 . 

Lately, the scalar field dark matter (SFDM) model has 
received much attention. When the scalar field contains a 
self-interaction this model is also called the Bose-Einstein 
Condensate (BEC) dark matter model, both names are used 
in the li terature and are interch angeable. The main idea is 
simple (|Guzman fc Matosl [ioool l. The SFDM model pro- 
poses that galaxies form by condensation of a scalar field 
(SF) with an ultra-light mass of the order of ~ 10 _22 eV. 
Therefore, when we mention the SFDM or BEC model in 
this paper we are describing a scalar field that conden- 
sates somehow and becomes Bose Einstein Condensate dark 
matter. From this mass it follows that the critical temper- 
ature of condensation T c ~ 1/m^/ 3 ~TeV is very high, 
thus, they form Bose-Einstein Condensates (BEC) drops 
very early in the universe. It has be en proposed that thes e 
drops are the haloes of galaxies fsee lMatos fc Urena 20011 ). 
i.e., that haloes are gigantic drops of SF. On the other side, 
big structures form like in the LCDM mode l, by hierarchy 
i|Matos fc Urena 200 J ; ISuarez fc Matosll201ll ), thus, all pre- 
dictions of the LCDM model at big scales are reproduced by 
SFDM. In other words, in the SFDM model the haloes of 
galaxies do not form hierarchically, they from at the same 
time and in the same way when the universe reaches the crit- 
ical temperature of condensation of the SF, in a similar way 
as water drops form in the clouds. From this it follows that 
all galaxies must be very similar because they formed in the 
same manner and at the same moment. Therefore, from this 
model we have to expect that there exist well formed galaxy 
haloes at bigger redshifts than in the LCDM model. In this 
model the scalar particles with this small mass are such that 
their wave properties avoid the cusp and reduce the high 
number of small satellites (|Hu.Barkana fc Gruzinovl |2000| ) 
which i s another problem that is still present i n the ACDM 
mode l (|see Robles fc Matos 2011 for a review! ; iKlvpin et al.l 
Il999l) .Summarizing, it is remarkable that with only one free 
parameter, the ultra- light scalar field mass (m^ ~ 10 _22 eV), 
the SFDM model fits: 

(i) The evolution of the cosmological densities 
l|Matos,Vazquez fc Maganall2009l ). 

(ii) The rotatio n curves of big galaxies ijHarkol [ioTTal ; 
iBernal et afll2008h and LSB galaxies. 



(iii) With this mass, the critical mass of collapse for a real 
scalar f ield is just 10 12 Mp, i. e., the one observed in galaxy 
haloes (|Alcubierre et al.ll2002l ). 

(iv) The scalar field has a natural cut off, thus the 
substructure in clusters of galaxies is avoided naturally. 



With a scalar field mass of m,j 



10 eV the amount 



of substructure is compatible with the observed one 
l|Matos fc Urena 200 ll ; ISuarez fc Matosll201ll ) . 

(v) We expect that SFDM forms galaxies earlier than the 
cold dark matter model, because they form BECs at a crit- 
ical temperature T c >> TeV. So if SFDM is right, we have 
to see big galaxies at big redshifts with similar features. 

(vi) And recently it has been demonstrated that SFDM 
haloes maintain satellite galaxies going around big galaxies 
for enough time to explain the existence of old stars in the 
satellites, provided that t he mass of the SF is just ~ 
10" 22 eV |Lora et alj|201ll ). 



The idea was fi rst considered by ISinl l|l994l ); 
IJi S. U. fc Sin S. J.I (Il994) and independently intro- 
duced bv lGuzman fc Matosl l|2000r i. In the BEC model, DM 
halos can be described in the non-relativistic regime, where 
DM halos can be seen as a Newtonian gas. If we consider a 
SF self-interaction, we need to add a quartic term to the SF 
potential, in this case the equ ation of state of the SF is that 
of a p olytope of index n=l (|Suarez fc Matosl 120111 ; lHarkd 
2011a). Different issues of BEC DM halos and the cosmo- 
logical behavior o f the BEC model have been studied in 
IColpi et all ll 19861) ; iGleiser et all dl988h: iBohmer fc Harko ' 



(2007); lHarkd (|2011al ); iMatos fc Guzman! (|200ll ); IChavani; 

hsm- 

It is a fact that any model trying to become a serious 
alternative to ACDM has to succeed not only in reproduc- 
ing observations in which the standard model fails but also 
has to keep the solid description at large scale. For this rea- 
son, our aim in this work is to test the BEC model with 
the t wo observations mention ed above, the cusp/core prob- 
lem l|Bohmer fc Harko 1120071 ) and the constant DM central 
surface density. In order to do this, we used the Thomas- 
Fermi approximation and a static BEC DM halo to fit ro- 
tation curves of a set of galaxies. However, so far there was 
no comparison between the density profile and the data, in 
this work we fill this blank by fitting rotation curves of 13 
high resolution low surface brightness (LSB) galaxies and 
additionally compare the fits to two characteristic density 
profiles 1) the cuspy Navarro-Frenk- White (NFW) profile 
that results from N-body simulations using ACDM and 2) 
the Pseudo Isothermal (PI) core profile. The comparison al- 
low us to show that the model is in general agreement with 
the data and with a core in the central region. For our sec- 
ond result we found that the meaning of a core is somewhat 
ambiguous. In order to clarify the meaning and unify the 
description, we propose a new definition for the core and 
core radius that allow us to decide when a density profile 
is cusp or core. Using this definition in the BEC model dis- 
cussed above we find that the BEC model can reproduced 
the constant value of fj,o and as a crosscheck we used the 
PI profile and find our results to be in very good agreement 
with observations. This argues in favor of the model and our 
definition. 

In section 2, we describe the density profiles we use to 
fit the rotation curves and we give our definition of the core 
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and core radius. In section 3, we fit the galaxy data and 
obtain the fitting parameters and the core radius for the PI 
and BEC density profiles. In section 4, we discuss our results 
and in section 5 we give our conclusions. 



2 DARK MATTER HALO DENSITY 
PROFILES 

In this section we provide the dark matter profiles that will 
be used in the analysis. In the last part of the section we 
briefly describe the usual meaning of the core and establish 
a new definition for it. 



2.1 BEC profile 

The case in which the dark matter is in the form of a 
static Bose-Einstein condensate and the number of DM par- 
ticles in the ground sta te is very large was considered in 
iBohmer fc Harko I (|2007f ). Following th is paper and assum - 
ing the Thomas- Fermi approximation (Dal fovo et al.|[l999h 
which neglects the anisotropic pressure terms that are rel- 
evant only in the boundary of the condensate, the system 
of equations describing the static BEC in a gravitational 
potential V is given by 



= -pVV 



(4) 



(5) 



(6) 



V 2 V = AivGp, 
with the following equation of state 
p(p) = U a p 2 , 

where Uq ~ ^\ ° , p is the mass density of the static BEC 
configuration and p is the pressure, as we are considering 
zero temperature p is not a thermal pressure but instead it 
is produced by the strong repulsive interaction between the 
ground state bosons. Assuming spherical symmetry and de- 
noting R as the radius at which the press ure and density are 
zero, the density profile takes the form ijBohmer fc Harkol 
l2007h 



Ps(r) 



b sin(fcr) 
Po 



kr 



(7) 



where k = y/Gm 3 /h 2 a = n/R and = p B (0) is the BEC 
central density, m is the mass of the DM particle and a is 
the scattering lenght. The mass at the radius r is given by 



m(r) 



sin(fcr 
kr 



cos(kr] 



(8) 



from here the tangential velocity Vb of a test particle at a 
distance r, is 



Vi(r) 



4nGpo ( sin(fer) 



cos(fcr" 



k 2 \ kr 

The logaritmic slope of a density profile is defined as 

_ djlog p) 
d(log r) 

using (7) in (10) it is obtained <|Harkdl2011bh 

a(r) ■■ 



nr i 


r nr \ 


l--cot( 









(9) 



(10) 



(11) 



Additio naly, the logar ithmic slope of the rotation curve is 
defined |Harkdl2011bT ) by 



from (9) we get 



1 - 



rfflog V) 
d(log r) 

ijir/R) 2 



1 - (nr/R) cot(nr/R) 



(12) 



(13) 



2.2 Pseudo Isothermal profile 

All the empirical core profiles that exist in the literature 
fit two parameters, a scale radius and a scale density. A 
characteristic profile of this type is 

= rTW' (14) 

this is the PI profile (|Begeman et al.lll99lh . Here R c is the 
scale radius and p^ 1 is the central density. The rotation 
curve is 



V(r)pi = ^AttGp^'R 2 ( 1 - — arctan 



(15) 



2.3 Navarro- Frenk- White profile 

The NFW profile emerges from numerical simulations 
that use only CDM and are based on the ACD M model 
ijDubinski et al.lll99ll : iNavarro etaDll996l . Il997l l. In addi- 
tion to this, we have chosen this profile because it is repre- 
sentative of what is called the cuspy behavior (a fa —1) in 
the center of galaxies du e to DM. The NFW density pro- 
file (INavarro et al.| [l997) and the rotation curve are given 
respectively by 



PNFw{r) 



Pi 



(r/Rs)(l + r/Rs) 2 



(16) 



V N fw{t) = ^/4-KGpiRh - 



/«( 1 + ^ 



r/Rs 



1 + r/R a 



(17) 

pi is related to the density of the universe at the moment 
the halo collapsed and R^ is a characteristic radius. 



2.4 Meaning of the core radius and cusp/core 
discrepancy 

In the large scale simulations that use collissionless cold dark 
matter the inner region of DM halos show a density distri- 
bution described by a power law p ~ r a with a » — 1, 
such behavour is what is now called a cusp. On the other 
hand, observations mainly in dwarf and LSB galaxies seem 
to prefer a central density going as p ~ r°. This discrep- 
ancy between observation and the CDM model receives the 
name of cusp/core problem. Among the empirical profiles 
most frequently used to describe th e constant density be - 
havior in these galaxies are the PI jBegeman et alJll99ll ). 
the iso thermal (lAthanassoula et al.|[l987 ) and the Burkert 
profile (|Burkertlll99'5f r Even though their behavior is similar 
in the central region and is specified by the central density 
fitting parameter, their second parameter called the core 
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radius does not represent the same idea. For instance, in 
the PI profile feg.(|14|)) we see that the core radius will be 
the distance in which the density is half the central den- 
sity. For the Burkert profile the core radius R b c urk will be 
whe n p burk (R b c url °) = po" rfc /4 and for an isothermal profile 
(I) (|Spano et all 120081 ') p 1 \R' C ) = p' /2 3/2 . Hence, we see an 
ambiguity in the meaning of the core radius, they get the 
same name but the interpretation depends on the profile. If 
we want to compare the central density of LSB galaxies with 
that of NFW, it usually suffice to have a qualitative com- 
parison, so far this is what we have been doing by fitting 
empirical profiles. However, high resolution rotation curves 
demand a more quantitative comparison. Indeed, if we want 
to test models by fitting RCs we have to know the specific 
meaning and size of the core, then we will be able to tell if 
a model is consistent with a cusp or not by making a direct 
comparison with the data. 

For this reason we ask ourselves the question: What and 
where is the core? To solve the previous ambiguity and to 
unify the concept for future comparisons, we found that a 
good definition for the core is a region where the density pro- 
file presents logarithmic slopes a — 1 and the core radius 
will be the radius at which the core begins, that is to say, 
for radius smaller than the core radius we will have a ^ — 1, 
this means that its value r' is determined by the equation 

a(r') = -1 (18) 

The advantages of this definition are that the interpreta- 
tion is independent of the profile chosen (also notice that 
it applies to the total density profile and is not restricted 
to that of DM) and in virtue of the same definition we can 
directly tell if a DM model profile is cored or cuspy. With 
our new definition the specific distance at which the core ra- 
dius occurs still depends on the profile chosen but now the 
physical interpretation is only one. In the following when we 
refer to both the core and core radius we adopt the previous 
interpretation. 

Applying the definition to ((7]) we get the core radius 
for the BEC profile Rb and for comparison we use (|14[) be- 
cause in turns out that the parameter R c corresponds to the 
core radius as defined above. Finally, fitting the NFW pro- 
file provides a direct comparison between a cusp and core 
and hence to the cusp/core problem. 



F568-3 F568-3 DM only 




02468 10 02468 10 

r (kpc) r (kpc) 



F583-1 F583-1 DM only 




02468 10 02468 10 

r (kpc) r (kpc) 



Figure 1. Contribution of the baryons to the rotation curve for 
F568-3 and F583-1. We denote observed data by black dots with 
error bars, dark matter with blue asterisks, the disk with cyan 
squares and the the gas with magenta squared boxes. We fit the 
figures on the left asumming the minimun disk hypothesis while 
for the ones on the right we substract in quadrature the baryons 
and fit only the dark matter. We notice that the barionic compo- 
nent is not dominant in the outer regions and that the difference 
in the fits is barely noticeable. 




02468 10 0123 
r (kpc) r (kpc) 

Figure 2. Logarithmic slope of the rotation curve (0) for three 
values of R. From left to right R = 6,8,10. The plots on the 
left show a common behavior always reaching zero before R. The 
figure on the right shows the region where a linear behavior in 
the RC speed is still valid, the (green) horizontal line bounds this 
region and the radius in which this boundary is reached was found 
to be r ~ 0.31-R in each curve. 



3 FITS AND DATA 

We see from © that the BEC model satisfies p ~ r° near the 
origin, but a priori this does not imply consistency with ob- 
served RCs. Therefore, we fit the profiles in section 2 to thir- 
teen high resolution observed RCs of a sample of LSB galax- 
ies. T he RCs were taken from a subsample of Ide Blok et alJ 
|200ll ). we chose galaxies that have at least 3 values within 
~ 1 kpc, not presenting bulbs and the quality i n the RC in 
Ha is good as defined in iMcGaugh S. S. et all (|200ll ). The 
RCs in this work omit galaxies presenting high asymmetries 
and included in the error bars are experimental errors in the 
velocity measurement, inclination and small asymmetries. 
Because the DM is the dominant mass component for these 
galaxies we adopt the minimum disk hypothesis which ne- 
glects baryon contribution to the observed RC. In order to 
show that in LSB and dwarf galaxies neglecting the effect of 



baryons is a good hypothesis, we include in Fig.l two repre- 
sentative examples (F568-3 and F583-1). For these galaxies 
we plot the contribution of the gas, disk and the dark matter 
separately. We did the fitting first considering the total con- 
tribution and then using only DM (marked as 14 and 15 in 
Table 1 and 2). We found no substantial difference in our val- 
ues, which can be seen from our results in Table 1 and 2. As 
the other galaxies belong to the cate gory of DM dominated 
galaxies as other authors have shown (|Kuzio de Narav et al.l 
120081 ; Ide Blok et al.ll200ll ). neglecting baryons in our analy- 
ses will not modify substantially our results. 

As the difference between a core and a cusp is most no- 
table only for data values inside 1 kpc and given that in the 
interval ~ 1 to 10 kpc the slopes of core and cusp profiles 
are very similar, which can lead to the wrong conclusion 
that cuspy halos are consistent with observations, we deter- 
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Figure 3. Observed LSB galaxy rotation curves and density profiles with the best halo fits. Below each RC is its density profile 
along with the fits. Shown are the PI (green dashed-line); the BEC (red solid line) and NFW (black dobble-dotted line) DM halo 
profiles, the observational data is drawn with error bars. The gray arrow denotes the best fit to the data within Ri and the vertical 
arrows denote the PI (blue) and BEC (magenta) core radius. 



mined the logarithm ic slope and the uncertainty following 
Ide Blok et aITl|200ll ) with the difference that we fit only the 
data within 1 kpc and that there is no need of an uncertain 
"break radius". 

In Table 1 we list the fitting parameters of the profiles 
of section 2, we also include the values of the logarithmic 
slope and its uncertainty, the value Ri denotes the nearest 
radius to lkpc where a data point is given. We obtain a by 
fitting values inside Ri, and we also report the core radius 
for the BEC profile Rb in order to compare it with R c . In 
Table 2 we report both the value of eq. ((3| for the BEC 
profile and the logarithm of eq. ([2]) for PI and BEC profiles. 

In Fig. 3, we show our fits to the RC data and the 
density profiles, also shown are the core radius in the BEC 
(magenta) and in PI (blue) profiles. The gray arrow is the 
fit that determines a, the size of the arrow denotes the fitted 
region and is bounded by 7?i. 



4 DISCUSSION 

In Fig. 2 we plot f3 using three different values of R. From 
this figure we see a common behavior. We have that f} is 
a decreasing function of r and is zero before R, which tell 
us that (9) always reaches a maximun before R. We also 
notice that there is a region in which Vb ~ r, we can take 
this region to be when 0.9 ^ /3 ^ 1. If we use /? = 0.9 as 



an upper bound for the region in which the linear behavior 
(Vb ~ r) remains valid, we obtain an upper bound radius 
of r « 0.31.R. This means that for values of r ^ 0.31i? we 
expect Vb ~ r. The latter can be used as a test to the 
BEC model by fitting the RCs and veryfing this solid-body 
behavior within the mentioned region. The fits of the RCs in 
Fig. 3 prove that the solid-body like behavior characterized 
by a linear increase of the velocity in the central region is 
satisfied by the BEC model, in fact, it is more consistent 
with the core PI and BEC profiles than the cuspy NFW. 

If we now turn to the density profiles, our fits within 
Ri give an average v alue of a=— . 27 ± 0.18 consistent 
with those obtained in Ide Blok et al.l (|2001 ) 
and with a= 



„, a=-0.2 ± 0.2 
-0.29 ± 0.07 reported by lOh Se-Heon et all 



l|201l[ ) analyzing 7 THINGS dwarf galaxies. The case of 
ESO1870510 might be considered to be consistent with 
NFW profile, however it is the innermost value that consid- 
erably decreases a, being an irregular galaxy more central 
data near the innermost region is required to discard the 
possibility of any violent event that might have caused such 
a slope value. 

The density profiles corresponding to the RCs fits are 
also shown in Fig. 3 for each galaxy. We see that the BEC 
fits slightly deviates for the farthest data points as a re- 
sult of the finite size of the radius R that is fixed by the 
same data. This discrepancy is due to the fact that the halo 
might be more extended than the value R. As a matter of 
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Table 1. Best-fit parameters of the profiles in section 2. 



(a) 



Label 


Galaxy 


Po 




„2 
Xpi 


Pa 


K 




1 J 

Kb 


„2 
Xb 






(M & pc °) 


(kpc) 




(Mqpc °) 


(kpc) 




(kpc) 




1 


ESO1200211 


0.0464 


0.57 


1.45 


0.0138 


2.92 




1.46 


0.44 


2 


ESO1870510 


0.0548 


0.96 


0.63 


0.0329 


2.93 




1.465 


0.2 


3 


ESO3020120 


0.0598 


1.89 


0.03 


0.0229 


8.95 




4.425 


0.92 


4 


ESO3050090 


0.0276 


2.09 


0.04 


0.0217 


4.81 




3.04 


0.1 


5 


ESO4880049 


0.1035 


1.62 


0.99 


0.0549 


5.36 




2.68 


0.79 


6 


U4115 


0.1514 


0.93 


0.10 


0.1438 


1.52 




0.76 


0.15 


7 


U11557 


0.0156 


5.37 


0.08 


0.0143 


9.38 




4.69 


0.07 


8 


U11611 


0.2065 


1.46 


0.14 


0.0771 


6.51 




3.225 


5.24 


9 


U11748 


1.678 


0.98 


1.43 


0.4205 


5.33 




2.665 


4.37 


10 


U11819 


0.0869 


3.03 


0.20 


0.0539 


8.86 




4.43 


1.09 


11 


U11583 


0.119 


0.64 


0.12 


0.0953 


1.43 




0.715 


0.84 


12 


F568-3 


0.0361 


3.01 


0.53 


0.0248 


8.78 




4.39 


0.11 


13 


F583-1 


0.0317 


2.6 


0.54 


0.019 


8.53 




4.26 


0.48 


14 


F568-3 (DM) 


0.0264 


3.44 


0.84 


0.0202 


8.96 




4.48 


0.25 


15 


F583-1 (DM) 


0.0279 


2.79 


0.58 


0.0177 


8.66 




4.33 


0.37 


(b) 


Label 


Galaxy 




Pi 


Rs 




a 


Aa 


Rl 






(xlO~ 3 M pc- 3 ) 


(kpc) 








(kpc) 


1 


ESO1200211 




2.45 


5.7 


0.24 


-0.04 





53 


0.95 


2 


ESO1870510 




0.761 


31.82 


0.05 


-1.09 





76 


1.13 


3 


ESO3020120 




2.65 


19.72 


0.32 


-0.2 





16 


1.4 


4 


ESO3050090 




0.0328 


705.67 


0.22 


-0.7 





03 


1.1 


5 


ESO4880049 




1.42 


52.27 


0.16 


-0.09 





06 


0.96 


6 


U4115 




0.139 


341.74 


0.78 


-0.27 





24 


0.89 


7 


U11557 




0.0108 


2849.65 


1.43 


0.2 





31 


0.9 


8 


U11611 




11.59 


14 


1.77 


-0.15 





31 


1.1 


9 


U11748 




204.58 


5.53 


3.41 


-0.38 





11 


1.1 


10 


U11819 




1.19 


101.24 


1.07 


-0.64 





13 


2.5 


11 


U11583 




0.136 


238.568 


0.81 


-0.2 





24 


0.95 


12 


F568-3 




0.378 


120.78 


2.36 


0.03 





02 


1.9 


13 


F583-1 




0.345 


102.349 


0.55 


-0.03 





04 


1.1 


14 


F568-3 (DM) 




0.0715 


515.68 


2.87 


0.34 





27 


1.9 


15 


F583-1 (DM) 




0.329 


102.349 


0.58 


0.02 





07 


1.1 



a fact, the more extended the "flat" outer region in the RCs 
the more conspicuous the discrepancy. The main reason of 
this comes from Fig. 2 where we infer that the rotation curve 
speed always presents a maximum value followed by a con- 
tirmos decrease, which means that to avoid the mentioned 
discrepancy we need that the BEC rotation curve profile 
remains approximately constant after its maximum. From 
Fig.l we see that the total rotation curve is dominated by 
the dark matter contribution, specially in the outer regions. 
Hence, unless the baryons become the dominant component 
in the outer regions, which does not seem to be observed, it 
is unlikely that adding the barionic contribution to the RCs 
in our galaxies will solve the discrepancy. 

Some solutions to keep the BEC rotation curve con- 
stant after its maximum include finite t emperature correc- 
tions to (fT4)l (|Harko fc Madarassvl 1201 il l, this suffice to al- 
leviate the latter problem in LSB galaxies and dwarfs but 
not for bigger galaxies. Other authors proposed including 



j pr 

vortex lattices (|Rmdler-Daller fc Shapircj201ll;IZinnerll201l h 
and adding more nodes (|Sinlll994l ; |ji S. U. fc Sin S. J.lll994r i 



in the sol ution of system Q and J5l). N evertheless, it can 
be shown (|Guzman fc Urena-L6pedl200&t ) that a systems of 
many nodes is unstable, therefore so far no final conclusion 
has been reached. 

When comparing the BEC and PI core radius we find 
a general difference of ~ 2 kpc, the core size in the PI pro- 
file is approximately 2 kpc smaller than the BEC core size, 
but the PI central density is larger. In U4115, U11557 and 
U 11583 both profiles are very similar which results in a sim- 
ilar core and central density values, this can also be taken 
as a consistency check for our core definition. 

Comparing the values of Rb in Table 1 we did not find 
a tendency to a common value. Assuming that the core ra- 
dius determines the transition where the DM distribution 
changes from the outer region to the inner constant central 
density, the lack of a unique value means that there is not a 
common radius at which this transition takes place. 

For our second test we use Rb to calculate • We have 
already seen that R c and Rb are generally different and Rb 
is not a fit parameter. Hence a priori Rb is not expected to 
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Table 2. Derived quantities from the parameters in Table 1. 



Label 


Galaxy 


log fM^± 




9 dm 
(xl0 _9 cms~ 2 ) 


I 




i AO _i_ n i q 
T.'+ii zll w. io 


i fin 4- n 


0.908 


2 




i 70 _|_ n 78 


i «q _|_ n 
l.UO zn u.oo 


2.15 


3 


FSO30201 20 


2 (It + ^2 


2 nn + o 51 

£i.\J\J _1_ \J . 'J X 


2.17 


4 


ESO3050090 


1.76 ± 0.59 


1.81 ± 0.60 


4.58 


5 


ESO4880049 


2.22 ± 0.11 


2.16 ± 1.16 


6.63 


6 


U4115 


2.14 ± 0.40 


2.03 ± 0.23 


4.93 


7 


U11557 


1.92 ± 0.31 


1.82 ± 0.54 


3.03 


8 


U11611 


2.97 ± 0.13 


2.40 ± 0.42 


11.3 


9 


U11748 


3.21 ± 1.98 


3.05 ± 0.78 


50.4 


10 


U11819 


2.42 ± 0.30 


2.37 ± 0.25 


10.8 


11 


U11583 


1.88 ± 0.25 


1.83 ± 0.53 


3.08 


12 


F568-3 


2.03 ± 0.35 


2.03 ± 0.53 


4.91 


13 


F583-1 


1.91 ± 0.23 


1.91 ± 0.65 


3.66 


14 


F568-3 (DM) 


1.95 ± 0.51 


1.95 ± 0.68 


4.09 


15 


F583-1 (DM) 


1.89 ± 0.13 


1.88 ± 0.59 


3.45 



a Both fiQ 1 and p^ units are Mq/pc 2 




5 






pi - 


_ 4 

_ O 
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cn 
o 
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01 23456789 10 11 12 13 
N 

Figure 4. Plot of log (^/Mopc~ 2 ) and log (/^j^ 7 /Mqpc~ 2 ) for 
each galaxy. N denotes the galaxy according to Table 1. Here we 
observe that these values remain approximately constant in both 
profiles, this also serves as a crosscheck for our definition of Rg 
in the BEC profile. The green dashed-line represents the mean 
values given in JT9} and ||20|I . 



correlate with p$ . However, with the values of Table 1 we 
obtain 



log(^/M pc- 2 ) = log p§ Rb = 2.05 ± 0.56 



log(/io /M Q pc 



logpo'Rc 



. ± 0.46 



(19) 
(20) 



for the average values in the BEC and PI profiles respec- 
tively. We see the excellent agreement of (I19|) with (|20l) that 



was used as a crosscheck and with {T} in which a much big- 
ger sample was used. The agreement has shown that the 
BEC model is capable of reproducing the constancy of the 
value p,o, something that because of the cuspy nature is not 
possible in the NFW profile. 

In Fig. 4 we plot the above values for each galaxy. We 
define the DM central surface density (mentioned in the in- 
troduction) for the BEC profile by 



< £ >o,dm= M<i 



(21) 



where M < r b is obtained from ((HJ evaluated at Rb- From 
Table 2 we see that for U11748 the value logfiQ~ PI is con- 
siderably above the rest and with the largest uncertainty. 
For this reason, in the following analysis we omit both, this 
value and the smallest one that corresponds to ESO1200211. 
Doing this we get an average value of eq. (19) given by 
dm>~ 191.35 Mqpc 2 , and for the acceleration 
felt by a test particle located in Rb due to DM only we have 
5dm(-Rs) ~ 5.2 xl0 -9 cms -2 broadly consistent with 1[3"]|. 

The fact that all galaxies present approximately the 
same order of magnitude in qdm(Rb) might suggest that 
Rb represents more than a transition towards a constant 
density, it can give us information about the close relation 
between DM and the baryons. Moreover, in view of the lack 
of a unique core radius, we can interpret the transition in 
DM distribution as an effect of crossing a certain accelera- 
tion scale instead of a radial length scale. Such interpretation 
reminds us that given in MOND but with the big difference 
that the acceleration scale found is for DM and is not a 
postulate of the model. 

To determine which interpretation causes the transition 
whether an acceleration scale or a length scale we will need 
to study the properties of larger samples of galaxies observed 
with the new telescopes. 



© 0000 RAS, MNRAS 000, 000-000 



Flat Central Density Profile and Constant DM Surface Density in Galaxies from Scalar Field Dark Matter 



5 CONCLUSIONS 

In this paper we find that the BEC model gives a constant 
density profile that is consistent with RCs of dark matter 
dominated galaxies. The profile is as good as one of the 
most frequently used empirical core profiles but with the 
advantage of coming from a solid theoretical frame. We fit 
data within 1 kpc and found a logarithmic slope a=— 0.27 
± 0.18 in perfect agreement with a core. It is important to 
notice that the cusp in the central regions is not a prediction 
that comes from first principles in the CDM model, it is a 
property that is derived by fitting simulations that use only 
DM. 

We established the ambiguity present in the usual inter- 
pretation of the core radius, we proposed a new definition 
for the core and core radius that takes away the ambigu- 
ity and that has a clear meaning that allows for a definite 
distinction when a density profile is core or cusp. 

Using our definition we find the core radius in the BEC 
profile to be in most cases over 2 kpc bigger than the core 
radius in the PI profile. We have assumed the DM particles 
are bosons and that a great number of them is in the ground 
state in the form of a condensate. This led us to good re- 
sults for our sample of galaxies, but it might be necessary 
to consider more than these simple hypotheses. 

As a second result and direct consequence of our core 
definition, we were able to obtain the constant value of no 
which is proportional to the central surface density. This 
result is one of several conflicts that jeopardize the current 
standard cosmological model. 

If we continue to observe even more galaxies with a core 
behavior, this model can be a good alternative to ACDM. 
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